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Reply to "Comment on 'The nature of slow dy- 
namics in a minimal model of frustration-limited 
domains'" by J. Schmalian, P.G. Wolynes and S. 

Wu. In a recent paperQ we reported simulation results 
for a minimal model of frustration-limited domains with 
energy functional 
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We concluded that in all but a very small unexplored 
portion of the phase diagram the Langevin equation cor- 
responding to Eq. ^ 
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with Gaussian random force r)(r,t), does not generate 
glassy dynamics. Although relaxation times do grow 
dramatically as a function of inverse temperature, we 
found no evidence for non-exponential multi-step relax- 
ation. A simple Hartree approximation is in fact suffi- 
cient to capture the computed dynamics almost quanti- 
tatively. The validity of a harmonic reference system, we 
reasoned, demonstrates that activated processes and the 
possibility of an entropy crisis are not essential for dy- 
namics at temperatures above (but not extremely close 
to) the fluctuation-induced first order transition temper- 
ature, r tr . We further noted that facile growth of long- 
range order upon cooling below r tr renders the study of a 
stationary, supercooled disordered state impossible. The 
relevance of our conclusions for the stripe glass scenario of 
Schmalian and Wolynes has recently been questioned by 
Schmalian et al. 2]- In particular, they point out that our 
simulations span a range of effective temperatures above 
that at which glassiness sets in, T c , as predicted by ei- 
ther the self-consistent screening approximation (SCSA) 
or dynamical mean- field theory (DMFT). 

First, we discuss the nature of the approximations 
used in and to calculate T c . Schmalian et al. 
describe SCSA and DMFT as more sophisticated ap- 
proaches than the mode-coupling approximation (MCA) 
analyzed in our work. Regardless of the level of so- 
phistication of these approximations, we chose to focus 
on the MCA because it is the clear analog of Gotze's 
mode coupling theory (MCT) for relaxation of density 
fluctuations in liquids 3]. Although MCT overestimates 
T c in model systems 0, simulations show a robust corre- 
lation between the predicted ergodic-nonergodic transi- 
tion and the onset of caging in supercooled liquids^ |(|. 
We found no hint of caging behavior in our simulations, 
even below the values of r for which the MCA predicts 
loss of ergodicity. Furthermore, the correspondence be- 
tween static correlations and dynamical relaxation in our 
minimal model is much stronger than that exhibited by 
supercooled liquids, and the non-Arrhenius temperature 




FIG. 1: Coarsening of a configuration whose structural corre- 
lations are consistent with the disordered solution of DMFT 
at r = —0.2, qo — 0.5, and A = 1. Field values 0(r) at the sur- 
face of our (periodically replicated) simulation cell are shown 
for times t = (top left), t = 500 (top middle), t = 1000 (top 
right), t = 2000 (bottom left), t = 5000 (bottom middle), and 
t = 10, 000 (bottom right). 



dependence of correlation times is coupled to exponen- 
tial, rather than multi-step, relaxation. These prominent 
features of the dynamics generated by (0), while interest- 
ing, are therefore not generically similar to those found 
in supercooled liquids. 

In the conclusion of our work we mentioned that 
growth of long-range lamellar order occurs even slightly 
below rtr, indicating a lack of stable glassy behavior 
at temperatures below those studied thoroughly in Q]. 
Here, we address the criticism of Schmalian et al. by 
presenting detailed results for r < r tr . For several dis- 
ordered configurations representative of the equilibrium 
state at t — —0.14, qo = 0.5 and A = 1, we have per- 
formed instantaneous quenches to r = —0.2, the tem- 
perature at which Schmalian et al. predict the onset 
of glassiness, and have followed the subsequent evolu- 
tion in time. Instead of the aging behavior typical of 
a supercooled liquid Q (weak time dependence of one- 
time quantities and relative stability of the disordered 
phase), coarsening occurs, i.e., the modulated order of 
the ground state rapidly appears in all regions of our 
simulation cell, and ordered domains gradually align. As 
a more explicit test of the stability of the glassy state de- 
scribed in 0, we have constructed initial configurations 
representative of the DMFT solution of Schmalian et al. 
for t — —0.2, drawing values for the Fourier components 
of 4>(r) from appropriate Gaussian distributions |8|. A 
typical trajectory originating at such a glassy state is 
depicted in Fig. 1. Again, the system steadily coarsens 
towards the ordered state. In Fig. 2 we show the peak of 
the structure factor, So(t), as a function of time for exam- 
ples of both types of initial conditions. After a short tran- 
sient, the two results merge and grow approximately with 
the power law expected from previous studies of models 
in the Brazovskii class^. These low-temperature results 
suggest that the topography of (JTJ in this region of config- 
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FIG. 2: Aging of the peak of the structure factor, So(t), for 
two trajectories at r = —0.2, go = 0.5, and X = 1. The dashed 
line corresponds to an initial configuration representative of 
equilibrium at r = —0.14. The solid line corresponds to an 
initial configuration consistent with the statistics of DMFT 
described in Q. Although these results are anecdotal, late 
stage growth is consistent with a power law So(t) cx t a with 
exponent a » 0.5. 

uration space differs from that envisioned by Schmalian, 
Wolynes and coworkers and that standard entropy cri- 
sis arguments must be modified in this context. In par- 
ticular, the aging dynamics we have computed contrast 
strongly with the aging of a supercooled liquid, imply- 
ing that the relevant region of the landscape does 
not closely resemble that of a supercooled liquid at low 
temperatures 0, E3 ■ 

We do not deny that experiments in high temperature 
superconducting materials reveal slow wandering of 
stripes. This behavior may indeed be consistent with 
coarsening in two dimensions, which permits very slow 
rearrangement of modulated structures. In addition 
disorder and other factors neglected in may also 
stabilize the quasi-disordered dynamics of hole-rich 
and hole-poor regions, perhaps giving rise to glassy 
dynamics rather than simple coarsening. One possibility 



is that more complicated terms should be included in 
||TJ that would stabilize the glassy state . Another 
possibility is that a stable glass state arises only in the 
limit of strong coupling (A — > oo). Investigation of these 
possibilities is certainly a worthy direction for future 
research. 
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